DEFENSE  DOCUMENTATION  CENTEr, 

FOR 

SCIENTIFIC  AND  TECHNICAL  INFORMATION 


CAMERON  STATION.  ALEXANDRIA,  VIRGINIA 


NOTICE:  When  government  or  other  drawings,  Bpeci- 
fications  or  other  data  are  used  for  any  purpose 
other  than  in  connection  with  a  definitely  related 
government  procurement  operation,  the  U.  S. 
Government  thereby  incurs  no  responsibility,  nor  any 
obligation  whatsoever;  and  the  fact  that  the  Govern¬ 
ment  may  have  formulated,  furnished,  or  in  any  way 
supplied  the  said  drawings,  specifications,  or  other 
data  is  not  to  be  regarded  by  implication  or  other¬ 
wise  as  in  any  manner  licensing  the  holder  or  any 
other  person  or  corporation,  or  conveying  any  rights 
or  permission  to  manufacture,  use  or  sell  any 
patented  invention  that  may  in  any  way  be  related 
thereto. 


BLANK  PAGES 
IN  THIS 
DOCUMENT 
WERE  NOT 
FILMED 


ABERDEEN  PROVING  GROUND,  MARYLAND 


DDC  AVAILABILITY  NOTICE 


Qualified  requesters  may  obtain  copies  of  this  report  from  DDC. 


The  findings  in  this  report  arc  not  to  be  construed 
as  an  official  Department  of  the  Army  position. 


BALLISTIC  RESEARCH  L  A  B  0  R  A  T  0  R  I  E 


REPORT  NO.  1220 


(Supersedes  BRIM  Report  No.  IU72) 
OCTOBER  1963 


APPROXIMATE  CIRCULAR  AND 
NON-CIRCULAR  OFFSET  PROBABILITIES  OF  HITTING 


Frank  E.  Grubbs 


Associate  Technical  Director 


RDT  &  E  Project  No.  1A013001A039 


ABERDEEN  PROVING  GROUND, 


MARY  j  A  .1  li 


BALLISTIC  RESEARCH  LABORATORIES 


REPORT  NO.  1220 


FEGrubbs/bJ 

Aberdeen  Proving  Ground,  Md. 
October  19&3 


APPROXIMATE  CIRCULAR  AND 
NON-CIRCUIAR  OFFSET  PROBABILITIES  OF  HITTING 

ABSTRACT 

For  equal  or  unequal  delivery  errors  and  an  offset  point  of  aim,  the 
chance  that  the  burst  point  of  a  warhead  occurs  within  a  given  distance  of 
a  selected  point  of  the  target  is  approximated  by  reference  to  the  Non-central 
Chi-square  distribution.  Offset  circular  and  non-circular  probabilities  of 
hitting  for  the  two  and  three  dimensional  cases  may  thus  be  approximated  with 
a  single,  straight-forward  and  rather  simple  technique  by  the  use  of  a  central 
Chi-square  distribution  with  fractional  number  of  degrees  of  freedom  or  a 
transformation  to  approximate  Normality.  Computations  of  probabilities  of 
hitting  are  illustrated  by  examples.  The  approximations  recommended  appear 
to  be  of  sufficient  accuracy  for  many  weapon  systems  evaluation  problems. 

This  report  supersedes  BRL  Memorandum  Report  No.  1472. 
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INTRODUCTORY  DISCUSSION 


A  problem  of  some  importance  to  the  weapon  systems  analyst  is  that  of 
finding  the  probability  of  hitting  a  circular  target  (two-dimensional  case) 
or  a  spherical  target  (three-dimensional  case)  whether  the  delivery  errors  ere 
equal  or  unequal  and  also  for  point  of  aim  or  center  of  impact  (C  of  I)  of  the 
rounds  either  coinciding  with  the  target  centroid  or  offset  from  it.  Moreover, 
in  practically  all  applied  problems  of  study,  the  analyst  does  not  require  great 
accuracy  in  such  probabilities.  Certainly,  two-decimal  accuracy  is  generally 
good  enough  end  in  most  cases  probabilities  of  hitting  in  error  by  as  much  as 
.02,  or  even  .03,  will  not  be  of  the  greatest  concern.  It  therefore  appears 
desirable  to  record  a  straight-forward,  unique  and  rather  simple  technique  for 
approximating  probabilities  of  hitting  for  all  of  the  various  cases  referred  to 
above. 

As  is  rather  well-known,  probabilities  of  hitting  of  the  type  discussed 
herein  are  intimately  connected  with  the  probability  distributions  of  quadratic 
forms  in  normal  variables,  and  a  great  amount  of  investigation  has  been  carried 
out  on  this  subject.  See  references  [12],  [9],  [5],  [15],  [13],  [6],  [14],  for 
example.  In  many  cases  for  published  theory  on  quadratic  forms,  however, 
computations  of  hit  probabilities  involve  infinite  series  expansions  or  special 
computer  programming  techniques,  matrix  notation  and  very  diverse  theoretical 
ramifications  which  the  average  weapon  systems  analyst  has  little  interest  in 
or  does  not  really  need  to  know  for  many  applied  problems.  Indeed,  some 
simplification  is  desirable. 

In  this  report  we  do  not  claim  originality  for  any  of  the  theory  used, 
rather  it  is  a  matter  of  Just  how  available  techniques  may  be  put  together  in 
order  to  obtain  a  single,  satisfactory  and  useful  method  for  approximating  the 
probabilities  desired. 


THE  PROBLEM  AND  RELATED  THEORY 

It  is  instructive  to  depict  the  problem  geometrically.  We  will  do  this 
for  the  three-dimensional  case  and  later  point  out  how  the  theory  may  be  extended 
to  either  the  two-dimensional  case  or  even  the  N-dimensional  case.  The  situation 
is  depicted  in  the  following  figure: 


In  the  figure,  the  target  center  (or  some  other  target  point  of  interest) 
is  located  at  (a^,  ct,,  a^),  the  center  of  impact  (C  of  I)  or  aim  point  of  the 
rounds  at  (p^,  pg,  p^)  and  the  actual  burst  position  of  the  warhead  at  (x^, 

Xg,  Xj).  We  assume  that  the  burst  position  coordinates  or  the  random  variables 
x^,  Xg,  Xj  are  independently  normally  distributed  with  mean  values  and  variances 
given  by 


E  (x±)  =  p^,  i  =  1,  2,  J 

(1) 

Var  (x.)  =  a2  -  a2.,  i  =  1,  2,  3 

(2) 

i 


o 


We  desire  the  probability  that  the  distance  J(x ^  -  o^)2  +  (xg  -  ag)2  +  (x,  -  a,)£ 


will  be  less  than  or  equal  to  a  chosen  value,  R.  Thus,  we  want 
Pr  [^  -  ax)2  +  (x2  -  dig)2  +  ( x 5  -  utj)2  <  R2} 


(3) 


which  implies  we  need  the  distribution  of  the  quadratic  form  on  the  left-hand 
side  of  the  Inequality  in  (3).  We  note  that  this  is  a  particular  quadratic 
form  not  involving  cross  product  teres  such  as  x^Xg,  etc. 

If  we  let 


2  2  2  2  2  2 
o  «  of  +  Og  +  Cj  and  v^  ■  ojo 


ui  "  (xi  "  ^i)/°i  -  ai)/°i 


(M 

(5) 


then  the  quadratic  form  of  (3)  may  be  rewritten  as 

3  3 


&  ,1<U1 '  *1>s  ■  &  • 01>i 


Hence,  we  have  the  equivalent  probability  statements 

3 


Pr 


[V*  <  R2/o2}  «  Pr{  ^  (xt  -  at)2  <  R2  j 


(6) 


(7) 


Now  the  u.  are  each  independently  normally  distributed  with  zero  mean  and 
^  2 

unit  variance.  Also,  the  (u^  +  a^)  are  each  distributed  as  non-central 

Chi-square  (Fisher,  1928)  with  one  degree  of  freedom  and  non-centrality 

2 

parameter  a^.  Generally,  if  we  define  the  non-central  X  as  Patnaik  [10], 


i.e. , 


,2 


2-- 


(8) 


then  the  probability  density  of  X  is  given  by  [10] 

.2 


P(X  ) 


X 

5 


X 

2 


v 


V  i*_L 

k  + 


,2'  2  f  1  ■  1  x1 


i)  221  i! 


(9) 
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where  \  *  £  a.  is  the  non- central  parameter  and  v  the  number  of  degrees 
i-l  2 

of  freedom.  We  have  the  ordinary  central  X  -distribution  when  all  the  aA  =  0. 

Hence ,  cur  generalized  quadratic  form  (6)  is  merely  the  sum  of  products 
2  2 

of  constants,  *  ay  a ,  and  non-central  Chi-squares,  i.e. ,  a  sum  of 
weighted  non-central  Chi-squares.  Incidentally,  it  is  obvious  that  (6) 
generalizes  to  H-dimensions,  although  for  weapon  systems  studies,  our  interest 
is  limited  to  the  two-  and  three-dimensional  coses,  of  course.  (For  the 
two-dimensional  case,  we  simply  use  1*1,2  instead  of  1*1,2, 3.) 


Patnaik  [10]  points  out  that  a  weighted  sum  of  non-central  Chi-squares  may 
be  approximated  by  fitting  its  first  two  moments  to  those  of  the  ordinary 
central  Chi-square,  X^.  Vow  the  mean  (m)  and  variance  (v)  of  ^  as  defined 
in  (5)  are  easily  found,  being 

3 

a  -  )  v.. 


(Wi  -  OL^f/a2 


(10) 


and 


(ii) 

Since  the  mean  and  variance  of  the  central  Chi-square  (X  )  are  n  and  2n 
degrees  of  freedom  (d.f.),  Respectively,  we  may  fit  the  sum  of  weighted 
non-central  Chi-squares  (X  )  to  the  central  X*"  by  noting  that  the  moments 
are  proper  when 


E 


and 


or  2m^ /v  ie  teken  equal  to  the  central  X“  and  the  equivalent  number  of  degrees 
2 

of  freedom  is  n  =  2m  /v. 


(%*)  -4 


(12) 


(13) 
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Thus,  in  sumnary  we  propose  to  take  the  quantity  X2  =  2m'|'2 /v  as  being 

p  O 

approximately  distributed  as  a  central  X  with  n  *  2m  /v  degrees  of  freedom 
and  m  and  v  are  functions  of  the  amount  of  offset  and  the  delivery  errors  as 
given  by  (10)  and  (ll).  The  computed  value  of  Chi-square  is  taken  as 
X2  =  2m^/v  =  2mR2/va2,  and  this  last  quantity  is  referred  to  a  table  of 
the  probability  Integral  of  the  Chi-square  distribution.  The  desired 

probability  is  P  *  P(X2/n)  *  J  f (X^JdX2.  (See,  for  example,  reference 

°  J  o 

HI. ) 

Alternatively,  the  desired  probability  is  given  also  by  the  following 
relation 

P  -  I  (X2/y£n,  n/2-1)  »  I  (R 2/o2Jv,  m2/v-l)  (14) 


where  I  (u,  p)  is  Karl  Pearson's  Incomplete  Gamma  Function  [16].  Moreover,  since 

2 

there  is  a  direct  relation  between  the  Gamna  distribution,  the  X  distribution, 
and  the  Pais son  distribution,  we  may  compute  the  desired  probability  from  a 

/  X* 

table  of  cumulative  Poisson  probabilities.  Letting  P  (c,  a)  be  the  chance  of 
c  or  more  (whole)  success  when  the  expected  number  of  successes  is  a,  then  the 
chance  we  seek  is  simply 


P  ■  P  (c,  a)  « 


(15) 


where  P  (c,  a) 


x 

r=c 


-a 


ir/r! 


(in  this  connection,  our  c  and  a  here  are  generally  fractional  values  so  that 
we  merely  interpolate  on  the  tabulated  whole  values  of  c  to  find  the  desired 
probability. ) 


This  notation  conforms  with  that  of  Molina's  Tables: 

Poisson's  Exponential  Binomial  Limit,  D.  Van  Mo&trand  Company,  Inc.,  New  York, 

WF. 
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Fijsiii-v,  : ri—i'-i*  vransformation  given  in  [5]  or  the  Wiison-dilferty 
2 

transformation  [17]  of  X  to  approximate  normal  variables  may  be  used.  The 
Fisher  transformation  is 


J2P  - 


J? n-1  ■  a  normal  variable  vith  zero  mean  and  unit 


variance,  or  for  our  case,  take 


(16) 

(17) 


as  a  normal  variable:  N  (0,  1),  then 


(18) 


or  ve  simply  refer  J UmR^/ a^v  - «Am2/v-l  to  a  table  of  the  cumulative, 
standardized  Normal  Integral. 

2 

The  Wilson-Hilf erty  transformation  [17]  of  X  to  a  normal  variable  is  more 
* 

accurate  and  generally  produces  rather  excellent  results  for  cases  of  practical 
interest.  The  Wilson-Hilf erty  transformation  is 


t  -  (l-2/9o) 

>/2/9n 


(19) 


where  t  is  approximately  normally  distributed  with  zero  mean  and  unit  variance: 
N  (0,  1).  Thus,  we  refer  the  quantity 


-  (1  -  v/9m2) 
Vv/9m2 


(20) 


to  a  table  of  the  cumulative  normal  integral  to  find  the  desired  probability. 


The  above  discussion  centered  around  the  trlvariate  case.  For  the 
bivariate  case,  it  is  obvious  that  we  delete  the  third  terms  of  (3)  and  (6) 
and  the  last  or  third  terms  in  each  of  the  summands  of  (10)  and  (11). 
Otherwise,  proceed  as  indicated. 


* 


Mathur  [7]  indicates  that  for  three  or  more  degrees  of  freedom,  the 

2 

Wilson-Hilf  erty  approximation  to  X  is  within  .007  in  probability. 
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In  the  above,  we  have  used  the  first  two  moments  (mean  and  variance)  of 
a  sum  of  weighted  Hon- central  Chi-sqycres  for  an  approximation  based  on  the 
central  Chi-square.  Pearson  [11]  points  out  that  a  three-moment  central 
Chi-square  approximation  for  the  distribution  of  Non-central  Chi-square  may 
be  quite  accurate.  Returning  to  the  quadratic  form  ft  in  Equation  (6),  we 
recall  that  its  mean  is  m  and  variance  is  v,  i.e.  Equations  (10)  and  (11).  The 
third  central  moment  of  ft  is  easily  found  to  be 
N 

h5  =  8  £  (vft  +  5 vx3  a12)  (21) 

i=l 


so  that  Pearson's  is  given  by 

.  N 


Oi  ■  “j2/*5  *  {8  I  (vi5  *  *i  ‘i2)}2/{2Z  ai2>}3 

i*=l  i=l 

(22) 

Thus,  following  Pearson's  suggestion  [11],  ire  would  propose  the  use  of  a  central 
Chi-square  approximation 


(X2,  -  ft  )/<Jzn  -  (ft  -  m)/yfv 


(23) 


where  the  number  of  decrees  of  freedom,  n  ,  which  is  i;oi. orally  fractional, 
would  be  obtained  from 


n  =  8/px 


(2M 


The  Wilson-Hilferty  transformation  to  approximate  normality  could  then  be 

applied  to  the  new  X  , ,  i.e.  we  take 

n 


X2,  =  (ft  -  m)  \Z2n*/v  +  n 


(25) 


and  refer 


t  = 


^,/n'  -  (l-2/9n’)j  /Jzj 9n' 


(26) 
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to  a  table  of  the  Normal  probability  integral.  (See  Example  3  below  for  a 
case  involving  Pearson's  three-moment  approximation.) 

DISCUSSION  ON  RELATED  WORK 

In  connection  with  the  work  of  this  paper,  Grad  and  Solomon  (5 ]  have 
investigated  very  thoroughly  the  exact  and  approximate  distributions  of 

k  2 

quadratic  forms  of  the  type  Q  =  £  ax  in  connection  with  probability  of 

i=l  1 

hitting  problems.  The  quadratic  form  herein  may,  of  course,  be  reduced  to 
k  2 

the  form  Q  =  £  a.x.  by  finding  the  appropriate  eigenvalues  involved. 
i=l  1  1 

Gilliland  [4 ]  ha3  investigated  series  expansions  for  the  integral  of  the 

non-circular  bivariate  normal  distribution  over  an  offset  circle.  Also, 

as  is  no  doubt  well-known,  Research  Memorandum  R-234  [6]  of  the  Rand 

Corporation  gives  "offset  circle"  probabilities  for  the  bivariate  circular 

case.  The  Rand  offset  circle  probabilities,  Q(r,p)  =  1  -  P(r,p),  give  for 

0^  =  cx  (i.e.  Og  =  c^)  one  minus  the  chance  that  a  sample  point  from  a  bivariate 

circular  Normal  distribution  will  lie  on  or  within  a  circle  of  radius  r,  with 

the  aim-point  offset  by  the  distance  p.  In  the  notation  of  this  paper,  a 

2  2  2 

very  good  approximation  of  P(r,p)  is  simply  Pr(X  <  mr  /vcr^)  where 
m  =  1  +  p2/2a2,  v  =1  +  p2/c2  ,  p2  =  (d^  -  u1)2  +  (d ^  -  u2)2  and  X2  =  Chi- 

2 

square  with  n  =  2m  /v  degrees  of  freedom.  We  have  checked  some  100 
representative  values  in  various  parts  of  the  Band  tables.  In  general  excellent 
agreement  was  obtained  and  the  greatest  discrepancy  for  the  approximation 
based  on  formula  (20)  was  .025  in  probability.  For  the  bivariate  case,  with 

a  -  a  ,  Helen  J.  Coon  of  the  Ballistic  Research  Laboratories  derived  the 

x  y 

formula  for  estimating  the  offset  circular  hit  probabilities  by  using  the 

o 

Wilson-Hilferty  transformation  of  the  non-central  X“. 

Di  Donato  and  Jarnagin  give  a  computer  program  for  numerical  evaluation 
and  tables  of  six  particular  probability  levels  (?  =  .OJ,  *20,  .JO,  .JO,  .90 
and  .95)  for  the  offset,  non-circular  bivariate  case.  A  check  indicates  good 
agreement  of  our  approximation  with  their  values.  (See  Example  3«) 
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Pachares  [9 1  gives,  for  the  distribution  of  a  definite  quadratic  form 
involving  central  normal  variables,  an  alternating  infinite  series  expansion 
which  converges  absolutely,  and  which  gives  an  upper  bound  for  the  cumulative 
distribution  function  if  one  stops  with  an  even  power  or  a  lower  bound  for 
an  odd  power  of  the  series  expansion.  Shah  and  Khatri  [15 ]  generalize  the 
result  of  Ifechares  to  the  case  of  non-central  normal  variables. 

Ruben  [13],  [14]  in  rather  elegant  fashion  studies  in  rigorous  detail  the 
probability  content  of  regions  under  spherical  normal  distributions,  the 
regions  including  half-spaces,  hyperspheres,  hypercones,  hjpercylinders, 
ellipsoids,  simplices  and  polyhedral  cones. 

Imhof  [6]  discusses  exact  and  approximate  methods  for  computing  the 
distribution  of  quadratic  forms  in  non-central  normal  variables.  We  believe 
that  our  recommended  computation  based  on  (20)  is  simpler  and  sufficiently 
accurate  for  weapon  systems  evaluations. 

Additional  Points  of  Interest 

222 

a)  If  the  target  is  subject  to  location  errors,  o_  ,  a-  ,  o„  in 

xl  x2  Jx3 
2  2  2 

addition  to  the  ever-present  round-to-round  variances  a  ,  a"  ,  a  about  the 

X1  x2  x* 

C  of  I,  we  may  put 


2 

°1 


2 

°2 


2  2 

°5  °T 
J  x3 


2  2  2  2 
c  ox  -t  4  0/j 

and  otherwise  proceed  as  uei'ore  in  i  lie  analysis. 

b)  With  the  theory  herein  it  is  easy  to  determine  the  approximate 
circular  or  spherical  probable  error,  even  for  unequal  delivery  errors. 

Equating  (20)  to  zero  and  solving  for  R,  we  obtain  the  radial  distance  which 
includes  50  percent  of  the  points,  i.e.,  the  circular  probable  error  (CEP)  or 
the  spherical  probable  error  (SEP).  For  example,  fur  the  bivariate  circular 
case  with  no  offset,  we  have  a1  =  a2  =  0/V2,  m=v-.l  and  we  find  R  =  32/27 

=  l.lfi5  as  compared  to  the  well-known  relation  CEP  =  *72  ln2  -  1.178  o. . 
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c)  It  is  noted  that  there  is  a  relationship  between  the  problem 
discussed  here  and  "coverage”  problems  in  weapon  systems  analysis.  The 
so-called  coverage  problem  involves  the  measure  of  two  -  and  three-dimensional 
random  sets  or  intersections  of  circles,  spheres,  etc.  for  various  delivery 
errors.  Indeed,  the  amount  or  fraction  of  coverage  for  one  circle  falling 

on  another  circle,  for  example,  is  determined  from  geometrical  considerations, 
whereas  the  probability  that  the  coverage  will  equal  or  exceed  a  given  value 
may  be  estimated  from  the  methods  discussed  herein. 

d)  Once  probabilities  of  hitting  are  available  as  computed  herein, 
lethality  or  vulnerability  data  may  be  easily  Included  in  order  to  obtain 
overall  "kill"  probabilities  for  a  weapon. 

c)  Finally,  it  is  remarked  that  in  many  weapon  assessment  cases  it  will 
not  be  of  great  concern  whether  one  uses  a  round  (or  spherical)  target  as  compared 
to  a  square  target,  since  available  vulnerability  data  or  lethality  data  or 
other  input  information  may  lead  to  some  lack  of  precision  anyway. 

Example  1 

Using  problem  1,  page  7^4,  of  Gilliland's  paper  [4]  for  a  bivariate 
example,  we  have 

R  =  0^,  =  C*2  ■  o,  ^  -  <J-|/5,  Mg  =  a^|/5»  Cg  *  Oj/ 2,  and 

2  2  2  „  2 
o  =  +  a2  ~  5^/4 

From  (10)  and  (11),  for  1=1,  2,  we  find  that  m  *  1.04  and  v  =  1.468.  Thu3, 

X2  =  2mR 2/a2v  =  2  (1.04)  (.8)/l.468  =  1.134  and  n  =  2m2/v  =  1.474.  From  the 

°  2 
Bioraetrika  Tables  [1],  wc  find  our  P  =  P  (XQ/n)  =  .570.  Equivalently,  using 

?  P  r~  2 

Karl  Pearson's  Incomplete  Comma  function,  I  (u,  p)  =  I  (R  /a  vv,  m  /v-l)  = 

p 

I  (.660,  -.263)  =  .575*  The  Poisson  approximation  gives  P  (n/2,  Xq/2)  = 

P  («737>  *567)  =  .558»  These  three  values  differ  no  doubt  because  of  the 
double  interpolation  required  in  the  various  tables. 
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For  the  Wilson- Hilferty  approximation  (19),  we  refer 


JjP/p-fllA1.-.  v/>2>  .  .1729 

Jv/9rfi 

to  a  table  of  the  cumulative  normal  integral  and  ve  get  P  «  .569,  which  is  of 

A 

sufficient  accuracy.  (The  Fisher  approximation  (lo)  to  X‘  gives  P  *  .544.) 

In  contrast  to  all  the  above  probabilities  Gilliland's  computed  value  of  P 

is  .577. 

Example  2 

If  we  consider  Gillixand's  problem  3,  page  766  of  [4],  the  data  are  the 
some  as  in  Example  1  except  now  ■  Cg,  and  we  find  m  ■  1.04  and  v  -  1.06. 

P  2  2  • 

Our  computed  value  of  Chi-square  is  *  2mR  /vo  -  .963’  with  an  equivalent 
number  of  degrees  of  freedom  n  ■  2m^/v  ■  2.003.  For  these  values,  we  find 
from  the  Biometrlka  Tables  that  P  =  P  (XQ/n)  =  .381.  Tables  of  the  Incomplete 
Gamma  Function  give  I  (u,  p)  *  I  (.481,  .00148)  -  .381  and  interpolation  in  the 
Poisson  Tables  gives  P  »  P  (n/2,  X^/2)  *  P  (1.002,  .481)  «  381.  The  Wllson- 
Hilferty  approximation  gives  P  *  .378  and  Gilliland's  computed  value  is  .382. 

Example  3 

From  the  Di  Donato  and  Jornagin  tables  (2  J ,  we  select  the  following  three 
cases. 

a)  For  a  probability  P  «  .50,  we  have  »  1,  Cg  *  3,  *  1,  Ug  =  ^ 

and  R  =  4.283.  Our  m  =  2.70,  v  =  7*44,  t  -•  -.0223  from  (20)  and  hence  P  =  .4921, 
which  is  in  error  by  only  -.00(39. 

b)  For  P  =  .90,  another  case  is  =  1,  *  6,  =  5,  |ig  =  ^0  and 

R  =  28.159*  Note  that  we  have  considerubie  offset  und  non-circularity.  In 
this  case,  our  m  =  12.486,  v  =  44.042,  t  =  1.291  from  (20)  and  P  =  .9017, 
so  that  the  error  is  only  +.0017. 

c)  For  a  very  extreme  case,  we  take  =  1,  Og  =  10,  ^  *  50>  Ug  = 

P  =  .20  and  R  =  49.696  from  the  Di  Donato  and  Jarnagin  tables.  We  have  m  =  25.762 
v  *  2.980,  t  =  -.7495  and  P  *  .2268,  so  that  the  error  is  +.0268. 
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Had  we  used  Pearson’s  three-rao^er.t  approximation  (25),  (24)  for  this  case, 
we  have  p,  -  .5044,  n  =  26.28,  X^t  =  20.78  and  t  =  -.7268  from  (26).  Thus, 

P  =  .2557  and  the  error  is  +.0337,  i.e.  somewhat  greater  and  at  the  expense  of 

more  involved  computation.  Nevertheless,  Imhof  [6]  indicates  that  Pearson's 

2  2 

three-moment  approximation  of  Central  X  to  weighted  non-central  X  may  be 
quite  accurate  generally  as  compared  to  this  single  case. 

Example  4 

Mr.  H.  L.  Merritt  of  the  Ballistic  Research  Laboratories  conducted  a 

Monte-Carlo  or  sampling  experiment  on  the  ORDVAC  electronic  computer  tc 

compare  probabilities  of  hitting  a  spherical  target,  based  on  200  trials 

for  each  probability,  with  thut  obtained  by  use  of  the  Wilson-Hilferty 

2 

transformation  (19)  of  X  to  approximate  normality.  The  computations  were 
programed  by  Mrs.  Jean  Jones  employing  the  FORAST  machine  language  or 
symbolic  coding. 

The  "target"  was  a  sphere  of  radius  unity,  located  at  the  origin. 

2  2  2  2  2  2  2 

Different  values  of  0^,  o2,  o^,  0  *  °i  +  °2  4  6114  various  C  of  I  locations 

were  employed  as  indicated  on  the  following  table.  The  sample  size  for  each 
machine  run  to  predict  probability  of  hitting  was  200  as  already  mentioned. 

The  sampling  results  are  indicated  in  the  ta!  le  below,  along  with  the  computed 
Wilson-Hilferty  values  of  probability. 
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Table  of  Commuted  and  Monte  Carlo  Probabilities 


C  of  I  Locations 

2 

0 

(0, 

0,  0) 

(1/2, 

0 
.  >• 

0 

(1,  0,  0) 

(2, 

0,  0) 

(3,  0, 

0) 

WH 

Ma 

WH 

Ma 

WH  Ma 

WH 

Ma 

WH 

Ma 

al  m  a2  m  °3 

.5 

.890 

.885 

.807 

.755 

.337  .353 

0 

0 

0 

0 

1.0 

.607 

.615 

•  530 

.505 

.261  .320 

.001 

.020 

0 

0 

2.0 

.313 

.320 

.307 

•  315 

.174  .205 

.004 

.035 

0 

0 

3.0 

.204 

.210 

.205 

.195 

.127  .135 

.006 

.040 

0 

.010 

°1  *  2o2  “  2ff- 

•  5 

.868 

.860 

.739 

.770 

.432  .435 

.019 

.040 

0 

0 

1.0 

.630 

.670 

.554 

.555 

.372  .405 

•  055 

.090 

.001 

0 

2.0 

.388 

.400 

.364 

.400 

.287  .290 

.092 

.125 

.011 

.035 

3.0 

.279 

.250 

.269 

.285 

.233  .235 

.104 

.105 

.023 

.050 

■ 

a? 

a 

ro 

Q 

VJI 

.5 

.878 

.890 

.738 

.735 

.378  .  380 

.005 

.015 

0 

0 

1.0 

.618 

.635 

.524 

.520 

.308  .355 

.021 

.040 

0 

0 

2.0 

.350 

.365 

.314* 

.393* 

.221  .265 

.044 

.050 

.002 

.010 

3.0 

.236 

.275 

.219 

.270 

.171  .210 

.053 

.070 

.006 

.010 

°1  "  2<t2  *  5ff3 

** 

•  5 

.740 

.774 

.739 

.750 

.447  .452 

.025 

.044 

0 

.001 

1.0 

.640 

.662 

.568 

.584 

.392  .408 

.068 

.087 

.002 

.009 

2.0 

.415 

.412 

.390 

.366 

.312  .306 

.109 

.122 

.015 

.029 

3.0 

.310 

.275 

.299 

.266 

.259  .223 

.123 

.120 

.030 

.046 

WH  =  Wllson-Hilf erty 

Ma  -  Machine  (ORDVAC)  ** 

#  For  this  block  1000  shots  for  each 

Indicates  greatest  discrepancy  probability  were  generated. 

TCie  results  are  generally  within  the  expected  sampling  error.  Since, 

* 

however,  the  greatest  discrepancy  occurred  for  the  values  marked  ,  then  10,000 
"shots"  were  generated  or  sampled  on  the  ORDVAC  for  this  condition  and  the 
experimental  probability  of  hitting  converged  to  .323  vs.  the  .314  for  the  WH 
approximation.  The  agreement  is  thus  very  good  indeed. 
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There  is  some  evidence  that  our  procedure  along  with  the  Wilson-Hilferty 
transformation  may  on  the  average  Just  slightly  underestimate  the  true 
probabilities  wc  seek. 


Example  5 

Olven  that  the  x.  or  shots  are  normally  distributed  with  zero  means  and 
2  1 

variances  ,  what  is  the  chance  that  a  point  (round)  will  lie  (or  hit)  on 
or  within  the  ellipsoid 


Now,  we  want  the  probability 


•  2  2  2 
where  u^  •  and  c^  »  /a  .  106611  6114  variance  of  0  ■  z  c^  u^ 

are  easily  found  to  be 


Hence,  we  may  take 

V 

o 

as  being  approximately  distributed  as  Chi-square  with  2m  /v  d.f. 
Using  the  Wilson-Hilferty  transformation,  we  would  refer 

t « -  ik  -  *Z2sli. 

*>/v/9m2 

to  a  table  of  the  standardized  normal  integral. 
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Example  6 

Suppose  the  x^  (burst  positions)  are  independently  normally  distributed 
with  unequal  means  (n^)  and  unequal  standard  deviations  (cm)  and  we  wish  to 
find  the  probability  of  hitting  on  or  within  the  offset  ellipsoid 

0  =  £  (*i  -  V2/ai2  5  1 

Now  putting  u^  =  (x^  -  and  k^  =  we  have 

0  =  ^  Vt  (ut  +  At)2  where  =  o^/a^ 

Moreover,  the  mean  and  variance  of  0  are  given  by 

m  =  E(0)  -  £  V.  (1  +  Aj2);  v  =  2  £  1*  (l  +  2k*) 

2  2 

and  we  may  take  2m 0/v  as  being  approximately  distributed  as  X  with  2m  /v 
degrees  of  freedom. 

Other  quadratic  forms  could  be  treated  similarly. 
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